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Abstract. The Mars Global Surveyor (MGS), Mars Express (MEX), and Venus Express (VEX) experienced several superior 
solar conjunctions. These conjunctions cause severe degradations of radio signals when the line of sight between the Earth and 
the spacecraft passes near to the solar corona region. The primary objective of this work is to deduce a solar corona model 
from the spacecraft navigation data acquired at the time of solar conjunctions and to estimate its average electron density. The 
corrected or improved data are then used to fit the dynamical modeling of the planet motions, called planetary ephemerides. 
We analyzed the radio science raw data of the MGS spacecraft using the orbit determination software GINS. The range bias, 
obtained from GINS and provided by ESA for MEX and VEX, are then used to derive the electron density profile. These 
profiles are obtained for different intervals of solar distances: from 12R Q to 215R Q for MGS, 6R Q to 152Rq for MEX, and 
form 12R Q to 154R Q for VEX. They are acquired for each spacecraft individually, for ingress and egress phases separately 
and both phases together, for different types of solar winds (fast, slow), and for solar activity phases (minimum, maximum). 
We compared our results with the previous estimations that were based on in situ measurements, and on solar type III radio 
and radio science studies made at different phases of solar activity and at different solar wind states. Our results are consistent 
with estimations obtained by these different methods. Moreover, fitting the planetary ephemerides including complementary 
data that were corrected for the solar corona perturbations, noticeably improves the extrapolation capability of the planetary 
ephemerides and the estimation of the asteroids masses. 
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1. Introduction 

The solar corona and the solar wind contain primarily ionized 
hydrogen ions, helium ions, and electrons. These ionized par- 
ticles are ejected radially from the Sun. The solar wind pa- 
rameters, the velocity, and the electron density are changing 
with time, radial distances (outward from Sun), and the so- 
lar cycles (Schwenn & Marsch 1990, 1991). The strongly tur- 
bulent and ionized gases within the corona severely degrade 
the radio wave signals that propagate between spacecraft and 
Earth tracking stations. These degradations cause a delay and 
a greater dispersion of the radio signals. The group and phase 
delays induced by the Sun activity are directly proportional to 
the total electron contents along the LOS and inversely propor- 
tional with the square of carrier radio wave frequency. 

By analyzing spacecraft radio waves facing a solar con- 
junction (when the Sun directly intercepts the radio signals 
between the spacecraft and the Earth), it is possible to study 



Table 1: Previous models based on in situ and radio science 
measurements (see text for detailed descriptions). 



Spacecraft 


Data type 


Author 


Mariner 6 and 7 


radio science 


Muhleman et al. (1977) 


Voyager 2 


radio science 


Anderson et al. (1987) 


Ulysses 


radio science 


Bird et al. (1994) 


Helios land 2 


in situ 


Bougeret et al. (1984) 


Ulysses 


in situ 


Issautier et al. (1998) 


Skylab 


Coronagraph 


Guhathakurta et al. (1996) 


Wind 


solar radio 


Leblanc et al. (1998) 




burst III 
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the electron content and to better understand the Sun. An 
accurate determination of the electron density profile in the 
solar corona and in the solar wind is indeed essential for 
understanding the energy transport in collision-less plasma, 
which is still an open question (Cranmer 2002). Nowadays, 
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mainly radio scintillation and white-light coronagraph mea- 
surements can provide an estimation of the electron density 
profile in the corona (Guhathakurta & Holzer 1994; Bird et al. 
1994; Guhathakurta et al. 1996, 1999; Woo & Habbal 1999). 
However, the solar wind acceleration and the corona heating 
take place between 1 to 10 R where in situ observations are 
not possible. Several density profiles of solar corona model 
based on different types of data are described in the literature 
(Table 1). The two viable methods that are generally used to de- 
rive these profiles are (Muhleman & Anderson 1981) (1) direct 
in situ measurements of the electron density, speed, and ener- 
gies of the electron and photons, and (2) an analysis of single- 
and dual-frequency time delay data acquired from interplane- 
tary spacecraft. 

We performed such estimations using MGS, MEX, and 
VEX navigation data obtained from 2002 to 201 1 . These space- 
craft experienced several superior solar conjunctions. This hap- 
pened for MGS in 2002 (solar activity maximum), for MEX in 
2006, 2008, 2010, and 2011 (solar activity minimum), and for 
VEX in 2006, 2008 (solar activity minimum). The influences 
of these conjunctions on a spacecraft orbit are severely noticed 
in the post-fit range and can be seen in the Doppler residuals 
obtained from the orbit determination software (see Figure 1). 

In section 2, we use Doppler- and range-tracking observa- 
tions to compute the MGS orbits. From these orbit determi- 
nations, we obtained range systematic effects induced by the 
planetary ephemeris uncertainties, which is also called range 
bias. For the MEX and VEX spacecraft, these range biases 
are provided by the ESA (Fienga et al. 2009). These range bi- 
ases are used in the planetary ephemerides to fit the dynamical 
modeling of the planet motions. For the three spacecraft, solar 
corona corrections were not applied in the computation of the 
spacecraft orbits. Neither the conjunction periods included in 
the computation of the planetary orbits. 

In section 3, we introduce the solar corona modeling and 
the fitting techniques that were applied to the range bias data. 
In section 4, the results are presented and discussed. In partic- 
ular, we discuss the new fitted parameters, the obtained aver- 
age electron density, and the comparisons with the estimations 
found in the literature. The impact of these results on plane- 
tary ephemerides and new estimates of the asteroid masses are 
also discussed in this section. The conclusions of this work are 
given in section 5. 

2. Data analysis of MGS, MEX, and VEX spacecraft 

2.1. Overview of the MGS mission 

The MGS began its Mars orbit insertion on 12 September 1997. 
After almost sixteen months of orbit insertion, the aerobrak- 
ing event converted the elliptical orbit into an almost circular 
two-hour polar orbit with an average altitude of 378 km. The 
MGS started its mapping phase in March 1999 and lost com- 
munication with the ground station on 2 November 2006. The 
radio science data collected by the DSN consist of one-way 
Doppler, 2/3 way ramped Doppler and two-way range obser- 
vations. The radio science instrument used for these data sets 
consists of an ultra-stable oscillator and the normal MGS trans- 



mitter and receiver. The oscillator provides the reference fre- 
quency for the radio science experiments and operates on the 
X-band 7164.624 MHz uplink and 8416.368 MHz downlink 
frequency. Detailed information of observables and reference 
frequency are given in Moyer (2003). 

2.1.1. MGS data analysis with GINS 

The radio science data used for MGS are available on the 
NASA PDS Geoscience website 1 . These observations were an- 
alyzed with the help of the GINS (Geodesie par Integrations 
Numeriques Simultanees) software provided by the CNES 
(Centre National d'Etudes Spatiales). GINS numerically inte- 
grates the equations of motion and the associated variational 
equations. It simultaneously retrieves the physical parameters 
of the force model using an iterative least-squares technique. 
Gravitational and non-gravitational forces acting on the space- 
craft are taken into account. The representation of the MGS 
spacecraft macro-model and the dynamic modeling of the orbit 
used in the GINS software are described in Marty et al. (2009). 

For the orbit computation, the simulation was performed 
by choosing two day data-arcs with two hours (approx. one 
orbital period of MGS) of overlapping period. From the over- 
lapping period, we were then able to estimate the quality of 
the spacecraft orbit determination by taking orbit overlap dif- 
ferences between two successive data-arcs. The least-squares 
fit was performed on the complete set of Doppler- and range- 
tracking data-arcs corresponding to the orbital phase of the mis- 
sion using the DE405 ephemeris (Standish 1998). To initialize 
the iteration, the initial position and velocity vectors of MGS 
were taken from the SPICE NAIF kernels 2 . 

The parameters that were estimated during the orbit fit- 
ting are (1) the initial position and velocity state vectors of the 
spacecraft, (2) the scale factors F^ and Fj for drag acceleration 
and solar radiation pressure, (3) the Doppler- and range resid- 
uals per data-arc, (4) the DSN station bias per data arc, and 
(5) the overall range bias per data-arc to account the geometric 
positions error between the Earth and the Mars. 

2.1 .2. Results obtained during the orbit computation 

In Figure 1, we plot the root mean square (rms) values of the 
Doppler- and range post-fit residuals estimated for each data- 
arc. These post-fit residuals represent the accuracy of the or- 
bit determination. To plot realistic points, we did not consider 
19% of the data-arcs during which the rms value of the post-fit 
Doppler residuals are above 15 mHz, the range residuals are 
above 7 m, and the drag coefficients and solar radiation pres- 
sures have unrealistic values. In Figure 1, the peaks and the 
gaps in the post-fit residuals correspond to solar conjunction 
periods. The average value of the post-fit Doppler- and two- 
way range residuals are less than 5mHz and lm, which ex- 
cludes the residuals at the time of solar conjunctions. 



1 http://geo.pds.nasa.gov/missions/ mgs/rsraw.html 

2 http://naif.jpl.nasa.gov/naif/ 
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Fig. 1: Doppler- and range rms of post-fit residuals of MGS for each two day data-arc. The residuals show the accuracy of the 
orbit determinations. The peaks and gaps in residuals correspond to solar conjunction periods of MGS. 



2.2. MEX and VEX data analysis 

The MEX and VEX radiometric data were analyzed done 
by the ESA navigation team. These data consist of two-way 
Doppler- and range measurements. These data sets were used 
for the orbit computations of MEX and VEX. However, de- 
spite their insignificant contribution to the accuracy of the orbit 
computation, range measurements are mainly used for the pur- 
pose of analyzing errors in the planetary ephemerides. These 
computations are performed with the DE405 (Standish 1998) 
ephemeris. The range bias obtained from these computations 
is provided by ESA, and we compared them to light-time de- 
lays computed with the planetary ephemerides (INPOP), ver- 
sion 10b, (Fiengaet al. 2011b), and the DE421 (Folkneret al. 
2008) ephemerides. For more details, see Fienga et al. (2009). 



2.2.1 . MEX: orbit and its accuracy 

Mars Express is the first ESA planetary mission for Mars, 
launched on 2 June 2003. It was inserted into Mars orbit on 
25 December 2003. The orbital period of MEX is roughly 6.72 
hours and the low polar orbit attitude ranges from 250 km (peri- 
center) to 1 1500 km (apocenter). 

The MEX orbit computations were made using 5-7-day 
track data-arcs with an overlapping period of two days between 
successive arcs. The differences in the range residuals com- 
puted from overlapping periods are less than 3 m, which rep- 
resents the accuracy of the orbit determination. As described 
in Fienga et al. (2009), there are some factors that have limited 
the MEX orbit determination accuracy, such as the imperfect 
calibrations of thrusting and the off-loading of the accumulated 
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Fig. 2: Systematic error (range bias) in the Earth-Mars and the Earth- Venus distances obtained from the INPOPlOb ephemeris: 
(Top panel) range bias corresponding to the MGS obtained for each two day data-arc; (middle and bottom panels) range bias 
corresponding to the MEX and the VEX. 



angular momentum of the reaction wheels and the inaccurate 
modeling of solar radiation pressure forces. 

2.2.2. VEX: orbit and its accuracy 

Venus Express is the first ESA planetary mission for Venus, 
launched on 9 November 2005. It was inserted into Venus orbit 
on 1 1 April 2006. The average orbital period of VEX is roughly 
24 hours and its highly elliptical polar orbit attitude ranges 
from 185km (pericenter) to 66500km (apocenter). However, 
data were almost never acquired during the descending leg of 
orbit, nor around periapsis (Fienga et al. 2009). 

The orbit was computed in the same manner as the MEX. 
The computed orbit accuracy is more degraded than MEX 



(Fienga et al. 2009). This can be explained by the unfavorable 
patterns of tracking data-arcs, the imperfect calibration of the 
wheel off-loadings, the inaccurate modeling of solar radiation 
pressure forces, and the characteristics of the orbit itself. The 
differences in the range residuals computed from overlapping 
periods are from a few meters to ten meters, even away from 
the solar conjunction periods. 



2.3. Solar conjunction: MGS, MEX, VEX 

The MGS, MEX, and VEX experienced several superior so- 
lar conjunctions. In June 2002, when the solar activity was 
maximum, the MGS SEP angle (see Figure 3) remained be- 
low 10° for two months and went at minimum to 3.325° ac- 
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cording to available data. For MEX, the SEP angle remained 
below 10° for two months and was at minimum three times: in 
October 2006, December 2008, and February 201 1. Similarly, 
the VEX SEP angle remained below 8° for two months and 
was at minimum in October 2006 and June 2008. The MEX 
and VEX superior conjunctions happened during solar min- 
ima. The influences of the solar plasma on radio signals during 
solar conjunction periods have been noticed through post-fit 
range and Doppler residuals, obtained during the orbit com- 
putations. Owing to insufficient modeling of the solar corona 
perturbations within the orbit determination software, no cor- 
rection was applied during the computations of the spacecraft 
orbit and range rate residuals. The peaks and gaps shown in 
Figure 2 demonstrate the effect of the solar conjunction on 
the range bias. The effect of the solar plasma during the MEX 
and VEX conjunctions on the radiometric data are described in 
Fienga et al. (2009), and for the MGS it is shown in Figure 1. 
The range bias (Figure 2) during solar conjunctions is used for 
deriving the electron density profiles of a solar corona model. 
These profiles are derived separately from the orbit determina- 
tion. 

3. Solar corona model 

3.1. Model profile 

As described in section 1, propagations of radio waves through 
the solar corona cause a travel-time delay between the Earth 
station and the spacecraft. These time delays can be modeled by 
integrating the entire ray path (Figure 3) from the Earth station 
(LEarth,/,,) to the spacecraft (L s / C ) at a given epoch. This model 
is defined as 



At = 



2cn 



1 r L -" 

x I 

:ri(/) JL Em „, s/n 



n cri (f) = 1.240 x 10 4 



N e (l) dL 



f 



(1) 



1 MHz 



where c is the speed of light, n c „ is the critical plasma den- 
sity for the radio carrier frequency /, and N e is an electron 
density in the unit of electrons per cm 3 and is expressed as 
(Birdetal. 1996) 



(2) 



where B and e are the real positive parameters to be de- 
termined from the data. R and I are the solar radius and radial 
distance in AU. F(9) is the heliolatitude dependency of the elec- 
tron density (Bird et al. 1996), where 6 represents the heliolat- 
itude location of a point along the LOS at a given epoch. The 
maximum contribution in the electron density occurs when / 
equals the MDLOS, p (see Figure 3 ), from the Sun. At a given 
epoch, MDLOS is estimated from the planetary and spacecraft 
ephemerides. The ratio of the MDLOS and the solar radii (R ) 
is given by r, which is also called the impact factor: 



P_ 



The electron density profile presented by Bird et al. (1996) 
is valid for MDLOS greater than 4R Q . Below this limit, turbu- 
lences and irregularities are very high and non-negligible. The 
solar plasma is therefore considered as inhomogeneous and ad- 
ditional terms (such r~ 6 and r~ 16 ) could be added to Equation 
2 (Muhleman et al. 1977; Bird et al. 1994). However, because 
of the very high uncertainties in the spacecraft orbit and range 
bias measurements within these inner regions, we did not in- 
clude these terms in our solar corona corrections. 

At a given epoch, the MGS, MEX, and VEX MDLOS al- 
ways remain in the ecliptic plane. The latitudinal variations in 
the coverage of data are hence negligible compared to the vari- 
ation in the MDLOS. These data sets are thus less appropriate 
for the analysis of the electron density as a function of heliolat- 
itude (Bird et al. 1996). Equation 2 can therefore be expressed 
as a function of the single-power-law (e) of radial distance only 
and be reduced to 



cm 



(3) 



On the other hand, at a given interval of the MDLOS in 
the ecliptic plane, Guhathakurta et al. (1996) and Leblanc et al. 
(1998) added one or more terms to Equation 3, that is 



4 (ir 



-d 



cm 



-3 



(4) 



with c ^ 4 and d = 2. 

To estimate the travel-time delay, we analytically integrated 
the LOS (Equation 1) from the Earth station to the spacecraft, 
using Equation 3 and Equation 4 individually. The analytical 
solutions of these integrations are given in Appendix A. 

In general, the parameters of the electron density profiles 
differ from one model to another. These parameters may vary 
with the type of observations, with the solar activity, or with the 
solar wind state. In contrast, the primary difference between the 
several models postulated for the electron density (r >4) is the 
parameter e (Equation 3), which can vary from 2.0 to 3.0. 

For example, the density profile parameters obtained by 
Muhleman et al. (1977) using the Man'«er-7radio science data 
for the range of the MDLOS from 5R to 100R o at the time of 
maximum solar cycle phase are 

(1.3 ± 0.9) x 10 8 (0.66 ± 0.53) x 10 6 _ 3 

Ne ~ ~~b + r 2.08±0.23 Cm " • 

The electron density profile derived by Leblanc et al. 
(1998) using the data obtained by the Wind radio and plasma 
wave investigation instrument, for the range of the MDLOS 
from about 1.3R© to 215R Q at the solar cycle minimum is 



N„ = 



0.8 x 10 s 0.41 x 10 7 0.33 x 10 6 



r 6 fA r 2 

Furthermore, based on in situ measurements, such as those 
obtained with the Helios 1 and 2 spacecraft, Bougeret et al. 
(1984) gave an electronic profile as a function of the MDLOS 
from 64.5R© to 215R Q as follows 



N, = 



6.14 

,2.10 



6 
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Similarly, Issautieret al. (1998) analyzed in situ measure- 
ments of the solar wind electron density as a function of helio- 
latitude during a solar minimum. The deduced electron density 
profile at high latitude (>40°) is given as 



,e ~ p2.003±0.015 

This Ulysses high-latitude data set is a representative sample 
of the stationary high-speed wind. This offered the opportu- 
nity to study the in situ solar wind structure during the minimal 
variations in the solar activity. As presented in Issautier et al. 
(1998), electronic profiles deduced from other observations in 
numerous studies were obtained either during different phases 
of solar activity (minimum or maximum), different solar wind 
states (fast or slow-wind), using data in the ecliptic plane (low 
latitudes). These conditions may introduce some bias in the es- 
timation of electronic profiles of density. 

Comparisons of the described profiles with the results ob- 
tained in this study are presented in section 4.2 and are plotted 
in Figure 9. 



Probe 




Fig. 3: Geometric relation between Earth-Sun-Probe. Where (3 
is the Earth-Sun-Probe (ESP) angle and a is the Sun-Earth- 
Probe (SEP) angle. 



3.2. Solar wind identification of LOS and data fitting 

As described in Schwenn (2006), the electronic profiles are 
very different in slow- and fast-wind regions. In slow-wind re- 
gions, one expects a higher electronic density close to the MNL 
of the solar corona magnetic field at low latitudes (You et al. 
2007). It is then necessary to identify if the region of the LOS is 
either affected by the slow-wind or by the fast- wind. To investi- 
gate that question, we computed the projection of the MDLOS 



on the Sun surface. We then located the MDLOS heliographic 
longitudes and latitudes with the maps of the solar corona mag- 
netic field as provided by the WSO 3 . This magnetic field is cal- 
culated from photospheric field observations with a potential 
field model 4 . 

However, zones of slow-wind are variable and not pre- 
cisely determined (Mancuso & Spangler 2000; Tokumaru et al. 
2010). As proposed in You et al. (2007, 2012), we took lim- 
its of the slow solar wind regions as a belt of 20° above and 
below the MNL during the solar minima. For the 2002 solar 
maximum, this hypothesis is not valid, the slow-wind region 
being wider than during solar minima. Tokumaru et al. (2010) 
showed the dominating role of the slow-wind for this entire pe- 
riod and for latitudes lower than ±70° degrees. 

An example of the MDLOS projection on the Sun's sur- 
face with the maps of solar corona magnetic field is shown in 
Figure 4. These magnetic field maps correspond to the mean 
epoch of the ingress and egress phases of the solar conjunc- 
tion at 20R Q . The dark solid line represents the MNL and the 
belt of the slow-wind region is presented by dashed lines. The 
two marked points give the projected locations of the MDLOS, 
ingress (a) and egress (t) at 2OR . For the entire MDLOS, this 
distribution of the slow- or fast-wind regions in the ingress and 
egress phases of solar conjunctions are shown in Figure 5. The 
black (white) bars in Figure 5 present the count of data sets 
distributed in the slow- (fast-) wind region. 

After separating the MDLOS into slow- and fast-wind re- 
gions as defined in Figure 5, the parameters of Equations 3 and 
4 are then calculated using least-squares techniques. These pa- 
rameters are obtained for various ranges of the MDLOS, from 
12R to 215R for MGS, 6R Q to 152?^ for MEX, and from 
12R to 154R for VEX. The adjustments were performed, for 
all available data acquired at the time of the solar conjunctions, 
for each spacecraft individually, and separately for fast- and 
slow-wind regions (see Table 2). 

To estimate the robustness of the electronic profile deter- 
minations, adjustments on ingress and egress phases were per- 
formed separately. The differences between these two estima- 
tions and the one obtained on the whole data set give the sen- 
sitivity of the profile fit to the distribution of the data, but also 
to the solar wind states. These differences are thus taken as the 
uncertainty in the estimations and are given as error bars in the 
Table 2. 

4. Results and discussions 

4.1. Estimated model parameters and electron density 

As described in section 3.2, we estimated the model parameters 
and the electron density separately for each conjunction of the 
MGS, MEX, and VEX. A summary of these results is presented 
in Table 2. The MDLOS in the unit of solar radii (R Q ) men- 
tioned in this table (column 5) represents the interval of avail- 
able data used for calculating the electronic profiles of den- 
sity. These profiles were then used for extrapolating the average 

3 http://wso.stanford.edu/ 

4 http://wso.stanford.edu/synsourcel.html 



Table 2: Solar corona model parameters and electron densities estimated from two different models using the MGS, MEX, and VEX range bias. The electron density out from 
the interval of the MDLOS is extrapolated from the given model parameters. 



S/C 



Year S.C 1 S.W.S 2 MDLOS 3 



N„ = Br 



N, =Ar 4 + Br 



B (x 10 6 ) 



Ne @ 20R a Ne @ 215R C 



B (xlO 6 ) A (x 10 8 ) Ne @ 20R e Ne @ 2i5R e 



MGS Aug 2002 Max 



S.W 12-215 



F.W 



0.51±0.06 2.00±0.01 1275+150 



1 1± 1.5 



0.52+0.04 



1300±100 



1 1 ± 1 



MEX Oct 2006 Min 



S.W 



F.W 



6-40 



1.90+0.50 2.54+0.07 942+189 



2.3+0.9 



0.30+0.10 0.16+0.02 850+263 



6+1 



MEX Dec 2008 Min 



S.W 



F.W 



6-71 



0.89+0.42 2.40+0.16 673+72 



2.3+1.3 



0.22+0.04 0.10+0.02 615+90 



5+1 



MEX Feb 20 11 Min 



S.W 40-152 0.52+0.10 2" 1300+25* 11 + 1 

F.W 6-60 1.70+0.10 2.44+0.01 1138+31 3.5+0.5 



0.52+0.10 1300+25 11 + 1 

0.33+0.02 0.24+0.04 975+60 7+1 



VEX Oct 2006 Min 



S.W 12-154 0.52+0.30 2.10+0.10 964+600 7+4 

F.W 12-130 1.35+1.10 2.33+0.30 1256+649 5+2 



0.40+0.28 
0.44+0.15 



1000+600 
1087+480 



9+6 
9+4 



VEX Jun 2008 Min 



S.W 12-154 1.70+1.50 2.50+0.50 950+625 3+2 

F.W 41-96 0.10+0.01 2" 250+25* 2+1 



0.31+0.20 
0.10+0.01 



775+450 
250+25* 



7+4 
2+1 



1 S.C: Solar cycle 

2 S.W.S: Solar wind state (F.W: fast-wind, S.W: slow-wind) 

3 



MDLOS: Minimum distance of the line of sight in the unit of solar radii (R Q ) 
Extrapolated value 
Fixed 
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Fig. 4: Solar corona magnetic field maps extracted from WSO at the mean epoch of 2OR ingress and egress. The dark solid line 
represents the magnetic neutral line. The dashed red lines correspond to the belt of the slow-wind region. The two marked points 
give the projected locations of the ingress (a) and egress (T) minimal distances at 20R Q . For the MGS 2002 solar conjunction, 
the hypothesis of a ±20° belt is not relevant (see section 4. 1 . 1) 



electron density at 215R (1AU). The period of the solar con- 
junctions, solar activities, and solar wind states are also given in 
columns 2, 3, and 4. Table 2 also contains the estimated param- 
eters of two different models: the first model from Bird et al. 
(1996) corresponds to Equation 3, whereas the second model is 
based on Guhathakurta et al. (1996) and Leblanc et al. (1998) 
and follows Equation 4 with c=4. Estimated model parameters 
for the slow- and fast-wind regions are presented in columns 6 
and 7. 



4.1 .1 . Mars superior conjunction 

The MGS experienced its superior conjunction in 2002 when 
the solar activity was maximum and the slow-wind region was 
spread at about +70° of heliolatitude (Tokumaruet al. 2010). 
Figure 6 represents the projection of the MGS MDLOS on the 
solar surface (black dots) superimposed with the 2002 synop- 
tic source surface map of solar wind speeds derived from STEL 
IPS observations extracted from Tokumaru et al. (2010). It sug- 
gests that, the MDLOS of the MGS exclusively remains in the 
slow-wind region. Respective estimates of the model parame- 
ters and of the electron densities are given in Table 2. From this 
table one can see that the estimates of the electron density from 
both models are very similar. Parameter e of Equation 3 is then 
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Fig. 5: Distribution of the latitudinal differences between MDLOS and MNL in the slow- and fast-wind regions during the ingress 
and egress phases of solar conjunctions. The black (white) bars present the slow- (fast-) wind regions as defined by +20° (>+20°) 
along the MNL during solar minima. 



estimated as 2.00+0.01, which represents a radially symmetri- 
cal behavior of the solar wind and hence validates the assump- 
tion of a spherically symmetrical behavior of the slow-wind 
during solar maxima (Guhathakurta et al. 1996). Whereas for 
Equation 4, the contribution of the r~ 4 term at large heliocen- 
tric distances (r > 12R ) is negligible compared to the r~ 2 term. 
Thus, the parameter A for these large heliocentric distances is 
fixed to zero and consequently gave similar results to Equation 
3, as shown in Table 2. 

In contrast, MEX experienced its superior conjunctions in 
2006, 2008, and 2011 during solar minima. The distribution 
of data during these conjunctions with respect to MDLOS are 
shown in Figure 7. From this figure, one can see that the MEX 
2006 (2008) conjunction corresponds to the fast- (slow-) wind 
region, whereas 2011 conjunction is a mixture of slow- and 
fast-winds. The estimated parameters of these conjunctions are 
given in Table 2. An example of the comparison between two 
models (Equation 3 and 4) is shown in Figure 8. This figure 
compares the electron density profiles obtained from the two 
models during the MEX 2008 conjunction. These profiles are 
extrapolated from 1R Q to 6R and from 71R Q to 215R Q - The 
upper triangles in Figure 8 indicate electronic densities ob- 
tained at 2OR and at the extrapolated distance of 215R Q (1 AU) 



with error bars obtained as described in section 3.2. As one sees 
in that figure, electronic profiles are quite similar over the com- 
putation interval till 71R and become significantly different 
after this limit. This suggests that the estimated parameters for 
both models are valid for the range of MDLOS given in Table 
2. Finally, for the MEX 201 1 conjunction, as shown in Figures 
5 and 7, the data are mainly distributed in the slow-wind (63%) 
during ingress phase and in the fast-wind (100%) during egress 
phase. Owing to the unavailability of slow-wind data near the 
Sun (MDLOS < 4OR ), we fixed e to 2 for Equation 3 (see 
Table 2). Moreover, the average electron density estimated for 
this conjunction is higher for the slow-wind than for the fast- 
wind and it is consistent with Tokumaruet al. (2010), which 
suggests that near the MNL, the electron content is higher than 
in the fast- wind regions. 

4.1 .2. Venus superior conjunction 

The VEX 2006 and 2008 conjunctions exhibit a mixture of 
slow- and fast-wind (Figure 7). These conjunctions occurred 
approximately at the same time as the MEX superior conjunc- 
tions. However, the limitations in the VEX orbit determination 
introduced bias in the estimation of the model parameters and 
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Fig. 6: 2002 synoptic source surface map of solar wind speeds derived from STEL IPS observations extracted from 
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the electron densities. This can be verified from the discrep- 
ancies presented in Table 2 for the VEX 2006 and 2008 con- 
junctions. Despite these high uncertainties, post-fit range bias 
corrected for the solar corona allows one to add complemen- 
tary data in the construction of the planetary ephemerides (see 
section 4.5). 



4.2. Comparison with other models 

Table 3 represents the estimated electron densities at 20R Q 
and 215R (1AU) from the various models. These models are 
representative of radio science measurements (Muhleman et al. 
1977; Anderson et al. 1987; Bird et al. 1994, 1996), in situ 
measurements (Bougeretet al. 1984; Issautieret al. 1997), and 
solar type III radio emission (Leblanc et al. 1998) measure- 
ments (Table 1). Table 3 and Figure 9 allow us to compare the 
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Authors Spacecraft Solar MDLOS Ne @ 20R o Ne @ 215Ro 



activity (el. cm 3) (el. cm 3) 



Leblanc et al. (1998) 


Wind 


Min 


1.3-215 


847 


7.2 


Bougeret et al. (1984) 


Helios 1 and 2 


Min/Max 


65-215 


890 


6.14 


Issautieret al. (1998) 


Ulysses 


Min 


327-497 


307* 


2.65+0.5* 


Muhleman et al. (1977) 


Mariner 6 and 7 


Max. 


5-100 


1231± 64 


9±3 


Birdetal. (1994) 


Ulysses 


Max 


5-42 


1700± 100 


4.7±0.415 


Anderson et al. (1987) 


Voyager 2 


Max 


10-88 


6650+ 850 


38±4 



* Mean electron density corresponds to latitude >40° 




MDLOS from 1-10 solar radii 




MDLOS from 10-215 solar radii 
Fig. 9: Comparison of different electron density profiles at different phases of the solar cycle from 1R to 215R (1AU). 



average electron density, obtained from different observations, 
made approximately during the same solar activity cycle. 

Figure 9 illustrates the comparisons of different electron 
density profiles, extrapolated from 1R to 215R Q . From this 
figure it can be seen that approximately all electron density 
profiles follow similar trends (oc r~ e , e varying from 2 to 3) 
until 1OR (panel E), whereas the dispersions in the profiles 
below 1OR (panel A) are due to the contribution of higher or- 
der terms, such as r~ 4 , r~ 6 or r~ 16 . 

In Table 3, we also provide the average electron density 
at 2OR and 215R , based on the corresponding models (if 
not given by the authors). The two individual electron den- 



sity profiles for ingress and egress phases have been given by 
Anderson et al. (1987) and Bird et al. (1994). To compare their 
estimates with ours, we took the mean values of both phases. 
Similarly, Muhleman et al. (1977) gave the mean electron den- 
sity at 215R estimated from round-trip propagation time de- 
lays of the Mariner 6 and 7 spacecraft. 

Table 3 shows a wide range of the average electron den- 
sities, estimated at 2OR and 215R during different phases 
of solar activity. Our estimates of the average electron density 
shown in Table 2 are very close to the previous estimates, es- 
pecially during solar minimum. The widest variations between 
our results and the earlier estimates were found during solar 
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conjunction using the model parameters given in Table 2. The 
error bars plotted in the figure correspond to the electron den- 
sity obtained at 2OR and 215R Q (see Table 2). 

Table 4: Statistics of the range bias before and after solar 
corona corrections. 



Pre-flt Post-fit 



s/c 


mean (m) 


cr(m) 


mean (m) 


cr(m) 


MGS, 2002 


6.02 


10.10 


-0.16 


2.89 


MEX, 2006 


42.03 


39.30 


0.85 


9.06 


MEX, 2008 


16.00 


20.35 


-0.10 


4.28 


MEX, 2011 


15.44 


19.20 


0.11 


6.48 


VEX, 2006 


5.47 


11.48 


-0.74 


6.72 


VEX, 2008 


3.48 


11.48 


-0.87 


7.97 



maxima and can be explained from the high variability of the 
solar corona during these periods. 

4.3. Post-fit residuals 

One of the objectives of this study is to minimize the effect 
of the solar corona on the range bias. These post-fit range bi- 
ases can then be used to improves the planetary ephemeris 
(INPOP). The pre-fit range bias represents the systematic er- 
ror in the planetary ephemerides during the solar conjunction 
periods. Figure 10 shows the pre-fit residuals over plotted with 
the simulated time delay (in units of distance), obtained from 
the solar corona model based on Equation 3. In contrast, the 
post-fit range bias represents the error in the ephemerides after 
correction for solar corona perturbations. 

From Figure 10 one can see that the systematic trend of the 
solar corona perturbations is almost removed from the range 
bias. The post-fit range bias of the VEX at low solar radii (es- 



pecially during the egress phase) is not as good as the MGS 
and MEX. This can be explained by the degraded quality of the 
VEX orbit determination (see section 2.2.2). The dispersion in 
the pre-fit and post-fit range bias is given in Table 4. The es- 
timated dispersions in the post-fit range bias are one order of 
magnitude lower than the dispersions in the pre-fit range bias. 
It shows a good agreement between the model estimates and 
the radiometric data. The corrected range bias (post-fit) is then 
used to improve the planetary ephemerides (see section 4.5). 



4.4. Model parameter dependency on the 
ephemerides. 

The range bias data are usually very important for constructing 
the planetary ephemerides (Folkner et al. 2008; Fienga et al. 
2011b). These measurements correspond to at least 57% of 
the total amount of data used for the INPOP construction 
(Fienga et al. 2009). Range bias data at the time of the solar 
conjunctions are not taken into account due to very high un- 
certainties (see Figure 2). Equations 3 and Figure 3 show the 
dependency of the density profile over geometric positions of 
the spacecraft (orbiting a planet) relative to the Earth and the 
Sun. The range bias used for this study includes the error in 
the geometric distance of Mars and Venus relative to the Earth. 
These errors are varying from one ephemeris to another and 
impact directly on the estimates of the mean electron density. 

Figure 11 illustrates the electron density estimated 
at 2OR for MGS, MEX and VEX using DE421 and 
INPOPlOb ephemerides. The dashed (DE421) and dotted- 
dashed (INPOPlOb) lines present the time limit up to which 
these ephemerides are fitted over range bias data. The curve 
represents the differences between INPOPlOb and DE421 esti- 
mations of Mars-Earth geometric distances. 

In particular, the Mars orbit is affected by the belt of aster- 
oids. The asteroid masses may cause a degradation in the esti- 
mates of the Mars orbit. Therefore, as one can see in Figure 11, 
the geometric differences between the INPOPlOb and DE421 
estimates of the Mars-Earth geometric distances are magni- 
fied from the extrapolation period onward. Hence, the electron 
densities estimated using DE421and INPOPlOb are consistent 
with each other within the error bars before the extrapolation 
period, whereas after this period, the DE421and INPOPlOb es- 
timates of the electron densities are quite different from each 
other. The sensitivity of the solar corona parameters and the 
electron densities, deduced from the analysis of the range bias, 
is then low as long as the computation is included in the time 
interval of the fit of the planetary ephemerides. However, out 
from the fitting time, the quality of these computations can 
be degraded by the extrapolation capability of the planetary 
ephemerides. Conversely, by fitting the planetary ephemerides 
(INPOP) including data corrected for the solar corona pertur- 
bations, some noticeable improvement can appear in the ex- 
trapolation capability of the planetary ephemerides and in the 
estimates of the asteroid masses (see section 4.5). 
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Fig. 10: Left panel: model-estimated solar corona over-plotted on the pre-fit residuals. Right panel: post-fit residuals after corona 
corrections. Top, middle, and bottom panels correspond to the MGS 2002, MEX 2008, and VEX 2006 conjunctions. 



4.5. Impact on planetary ephemerides. 

As one can see in Figure 10, correct the effects induced 
by the solar corona on the observed Mars-Earth distances 
is significant over some specific periods of time (during so- 
lar conjunctions). We aim to estimate the impact of this im- 
portant but time-limited improvement of the measurements 
of interplanetary distances on the construction of the plane- 
tary ephemerides. To evaluate any possible improvement, we 
produced two ephemerides, INPOPlOc and INPOPlOd, both 
fitted over the same data set as was used for the construc- 
tion of INPOPlOb (Fiengaetal. 2011b). This data set con- 
tains all planetary observations commonly used for INPOP 
(see Fiengaet al. (2009, 2011a)), including the MGS data ob- 
tained in section 2.1.1 and the MEX and VEX range bias pro- 
vided by ESA. These newly built ephemerides are based on the 



same dynamical modeling as described in Fienga et al. (2009, 
2011a). However, INPOPlOc is estimated without any solar 
corona corrections on the MGS, MEX and VEX range bias, 
and INPOPlOd includes the solar corona corrections evalu- 
ated in the previous sections. The selection of the fitted aster- 
oid masses and the adjustment method (bounded value least- 
squares associated with a priori sigmas) are the same for the 
two cases. The weighting schema are also identical. The differ- 
ences remain in the quality and the quantity of the range bias 
used for the fit (one corrected for solar plasma and one not) and 
in the procedure selecting the data actually used in the fit. 

For INPOPlOc, about 119901 observations were selected. 
Of these, 57% are MGS, MEX, and VEX range bias data that 
are not corrected for solar corona effects. Based on a very 
conservative procedure, observations obtained two months be- 
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Fig. 11: Variation of the average electron density at 20R Q and 215R using the DE421 and INPOPlOb ephemeris. The dotted- 
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fore and after the conjunctions were removed from the fitted 
data sample. This strategy leads to removal of about 7% of 
the whole data set, which represents 14% of the MGS, MEX, 
and VEX observations. For INPOPlOd, thanks to the solar 
corona corrections, only observations of SEP smaller than 1 .8 
degrees were removed from the data sample. This represents 
less than 1 % of the whole data sample. The estimated accuracy 
of the measurements corrected for the solar plasma is 2.4 me- 
ters when observations not affected by the solar conjunctions 
have an accuracy of about 1.7 meters. By keeping more obser- 
vations during solar conjunction intervals, the number of data 
with a good accuracy is then significantly increased. 

Adjustments of planet initial conditions, mass of the sun, 
sun oblateness, mass of an asteroid ring, and the masses of 289 
asteroids were then performed in the same fitting conditions as 
INPOPlOb. 

No significant differences were noted for the evaluated pa- 
rameters except for the asteroid masses. 

For the masses estimated both in INPOP 10c and 
INPOPlOd, 20% induce perturbations bigger than 5 meters on 
the Earth-Mars distances during the observation period. The 
masses of these 59 objects are presented in Tables 5, 6, and 
7. Within the 1-cr uncertainties deduced from the fit, we no- 
tice 10 (17%) significant differences in masses obtained with 



INPOPlOc and INPOPlOd, quoted with a "*" in column 5 of 
the two tables, 7 (12%) new estimates made with INPOPlOd, 
noted N in column 5, and 6 (10%) masses put to in INPOPlOd 
when estimated in INPOPlOc, marked with in the fifth col- 
umn. 



Table 8 lists masses found in the literature compared with 
the quoted values of Tables 5, 6, and 7. In this table, 80% of 
the INPOPlOd estimates have a better consistency with the 
values obtained by close encounters than the one obtained 
with INPOPlOc. Of these, the new estimates obtained with 
INPOPlOd for (20), (139) and (27) agree well with the values 
found in the literature. For (45) Eugenia, the INPOPlOc value 
is closer to the mass deduced from the motion of its satellite 
(Marchis et al. 2008) even if the INPOPlOd value is still com- 
patible at 2-cr. For (130) Elektra, the INPOPlOc and INPOPlOd 
estimated values are certainly under evaluated. Finally, one can 
note the systematic bigger uncertainties of the INPOPlOd es- 
timates. The supplementary data sample collected during the 
solar conjunctions that has 30% more noise than the data col- 
lected beyond the conjunction can explain the degradation of 
the uncertainties for the INPOPlOd determinations compared 
to INPOPlOc. 
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By correcting the range bias for the solar corona effects, we 
added more informations related to the perturbations induced 
by the asteroids during the conjunction intervals. 

In principal, during the least-squares estimation of the as- 
teroid masses, the general trend of the gravitational perturba- 
tion induced by the asteroid on the planet orbits should be de- 
scribed the most completely by the observable (the Earth-Mars 
distances) without any lack of information. In particular, for 
an optimized estimation, the data sets used for the fit should 
include local maxima of the perturbation. 

However, it could happen that some of the local maxima 
occur during the solar conjunction intervals. One can then ex- 
pect a degradation of the least-squares estimation of the per- 
turber mass if no solar corrections are applied or if these in- 
tervals are not taken into account during the fit. To estimate 
which mass determination can be more degraded than another 
by this window effect, we estimated L, the percentage of local 
maxima rejected from the INPOP 1 Oc fit in comparison with the 
INPOPlOd adjustment including all data sets corrected for so- 
lar plasma. L will give the loss of information induced by the 
rejection of the solar conjunction intervals in terms of highest 
perturbations. 

The L criteria are given in column 7 of Tables 5, 6, and 
7. As an example, for (24) Themis one notes in Table 8 the 
good agreement between the close encounter estimates and the 
INPOPlOd mass determination compared with INPOPlOc. On 
the other hand, based on the L criteria, 36% of the local max- 
ima happen near solar conjunctions. By neglecting the solar 
conjunction intervals, more than a third of the biggest perturba- 
tions are missing in the adjustment. This can explain the more 
realistic INPOPlOd estimates compared with INPOPlOc. 

We also indicate in Tables 5, 6, and 7 if important con- 
straints were added in the fit (column 8). In these cases, even 
if new observations are added to the fit (during the solar con- 
junction periods), there is a high probability to obtain a sta- 
ble estimates of the constrained masses as for the biggest per- 
turbers of Table 5. For the other mass determinations, one can 
note a consistency between high values of the L criteria and 
the non-negligible mass differences between INPOPlOc and 
INPOPlOd. 




Fig. 12: MEX extrapolated residuals estimated with INPOPlOd 
(light dots) and INPOPlOc (dark dots). 



By improving the range bias residuals during the solar con- 
junction periods, we then slightly improved the asteroid mass 
determinations. 

Estimates of residuals for data samples not used in the 
INPOP fit and dated after or before the end of the fitting interval 
are currently made to evaluate the real accuracy of the planetary 
ephemerides (Fienga et al. 201 la, 2009). To estimate if the use 
of the solar corona corrections induces a global improvement 
of the planetary ephemerides, the MEX extrapolated residuals 
were computed with INPOPlOc and INPOPlOd. As one can 
see in Figure 12, the INPOPlOd MEX extrapolated residuals 
show a better long-term behavior compared with INPOPlOc 
with 30% less degraded residuals after two years of extrapola- 
tion. 

Supplementary data of the MEX and VEX obtained dur- 
ing the first six months of 2012 would confirm the long-term 
evolutions of the INPOPlOd, INPOPlOc and INPOP 10b. 

This improvement can be explained by the more realistic 
adjustment of the ephemerides with denser data sets (7%) and 
more consistent asteroid mass fitting. 

5. Conclusion 

We analyzed the navigation data of the MGS, MEX, and VEX 
spacecraft acquired during solar conjunction periods. We esti- 
mated new characteristics of solar corona models and electron 
densities at different phases of solar activity (maximum and 
minimum) and at different solar wind states (slow and fast). 
Good agreement was found between the solar corona model 
estimates and the radiometric data. We compared our estimates 
of electron densities with earlier results obtained with differ- 
ent methods. These estimates were found to be consistent dur- 
ing the same solar activities. During solar minima, the elec- 
tron densities obtained by in situ measurements and solar radio 
burst III are within the error bars of the MEX and VEX es- 
timates. However, during the solar maxima, electron densities 
obtained with different methods or different spacecraft show 
weaker consistencies. These discrepancies need to be investi- 
gated in more detail, which requires a deeper analysis of data 
acquired at the time of solar maxima. 

The MGS, MEX, and VEX solar conjunctions data allows 
us to analyze the large-scale structure of the corona electron 
density. These analyses provide individual electron density pro- 
files for slow- and fast-wind regions during solar maxima and 
minima activities. 

In the future, planetary missions such as MESSENGER 
will also provide an opportunity to analyze the radio-science 
data, especially at the time of maximum solar cycle. 

We tested the variability caused by the planetary 
ephemerides on the electron density parameters deduced from 
the analysis of the range bias. This variability is smaller than 
the 1-cr uncertainties of the time-fitting interval of the plan- 
etary ephemerides but becomes wider beyond this interval. 
Furthermore, data corrected for solar corona perturbation s were 
used for constructing of the INPOP ephemerides. Thanks to 
these supplementary data, an improvement in the estimation of 
the asteroid masses and a better behavior of the ephemerides 
were achieved. 
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Table 5: Masses of the 59 asteroids that induce perturbations greater than 5 meters on the Earth-Mars distances during the period 
of observations (Kuchynka et al. 2010). Columns 2 and 3 give the values of the masses with the 1-cr uncertainties obtained with 
INPOPlOc and INPOPlOd. The differences between INPOPlOc and INPOPlOd values are given in column 4. Column 5 indicates 
if the INPOPlOd masses are newly determined (quoted N), set equal to when estimated by INPOPlOc (quoted 0), significantly 
different from INPOPlOc determinations (*). The maximum impacts of the perturbations are given in column 5. In column 7 is 
given the L criteria, the rate of local maxima rejected from the INPOPlOc fit but included in the INPOPlOd. In the last column 
we provide the status of the imposed constraints: S for strong and N for normal. 



IAU designation 


INPOPlOc 


INPOPlOd 


Diff 


S 


Impact 


L 




number 


10 12 xM G 


10 12 xM Q 


10 12 x M Q 




m 


% 




4 


130.109 ± 0.716 


130.109 ±0.983 


0.000 




1198.953 


20.5 


S 


1 


467.267 ± 2.047 


467.267 ± 2.437 


0.000 




793.741 


17.3 


S 


2 


103.843 ± 1.689 


102.654 ± 1.933 


1.189 




146.270 


11.8 


s 


324 


5.723 ±0.531 


5.723 ±0.611 


0.000 




93.536 


1.0 


s 


10 


43.513 ±3.300 


43.513 ± 3.877 


0.000 




77.003 


15.9 


s 


19 


3.884 ± 0.447 


3.450 ± 0.526 


0.435 




59.069 


13.8 


N 


3 


11.793 ±0.714 


11.793 ±0.803 


0.000 




55.639 


0.6 


s 


704 


19.217 ±2.315 


19.217 ± 2.869 


0.000 




34.492 


7.4 


s 


532 


2.895 ± 1.043 


2.895 ± 1.093 


0.000 




32.714 


2.3 


s 


9 


3.864 ±0.613 


3.063 ± 0.665 


0.801 




29.606 


20.6 


N 


7 


5.671 ±0.512 


5.367 ±0.591 


0.305 




27.822 


13.9 


s 


29 


7.629 ± 1.067 


7.227 ± 1.225 


0.402 




26.673 


2.9 


s 


24 


7.641 ± 1.596 


2.194 ± 1.775 


5.447 


* 


26.131 


36.0 


N 


31 


3.256 ± 2.034 


4.411 ±2.050 


1.155 




23.466 


24.1 


s 


15 


13.576 ±0.939 


13.576 ± 1.264 


0.000 




21.555 


20.6 


s 


6 


7.084 ± 0.822 


7.084 ± 1.048 


0.000 




21.150 


7.4 


s 


11 


3.771 ±0.976 


3.771 ± 1.110 


0.000 




17.301 


31.9 


s 


139 


0.000 ± 0.000 


3.579 ± 0.595 


3.579 


N 


16.687 


32.0 


N 


747 


4.129 ±0.841 


6.805 ± 1.089 


2.676 


* 


15.937 


31.6 


N 


105 


3.111 ±0.556 


3.111 ±0.745 


0.000 




15.196 


4.5 


N 


20 


0.000 ± 0.000 


1.921 ±0.661 


1.921 


N 


14.763 


39.9 


N 


372 


12.365 ± 2.676 


12.365 ± 2.990 


0.000 




13.796 


19.0 


S 
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Table 6: Same as Table 5 



IAU designation 


INPOPlOc 


INPOPlOd 


Diff 


S 


Impact 


L 




number 


10 12 x M Q 


10 12 x M G 


10 12 x M Q 




m 


% 




8 


3.165 ± 0.353 


3.325 ± 0.365 


0.159 




12.664 


17.7 


S 


45 


3.523 ± 0.819 


1.518 ± 0.962 


2.005 


* 


11.790 


21.0 


N 


41 


3.836 ± 0.721 


2.773 ± 0.977 


1.063 




11.568 


15.2 


N 


405 


0.005 ± 0.003 


0.006 ± 0.003 


0.001 




11.378 


21.2 


N 


511 


9.125 ± 2.796 


9.125 ± 3.138 


0.000 




10.248 


20.5 


S 


52 


8.990 ± 2.781 


8.990 ± 3.231 


0.000 




9.841 


3.0 


S 


16 


12.613 ± 2.286 


12.613 ± 2.746 


0.000 




9.701 


8.9 


s 


419 


1.185 ± 0.461 


0.425 ± 0.398 


0.760 




9.585 


10.2 


N 


78 


0.026 ± 0.016 


0.024 ± 0.016 


0.002 




9.389 


9.8 


N 


259 


0.092 ± 0.002 


0.006 ± 0.003 


0.086 




9.222 


31.4 


N 


27 


0.000 ± 0.000 


1.511 ± 0.982 


1.511 


N 


9.146 


29.5 


N 


23 


0.000 ± 0.000 


0.093 ±0.156 


0.093 


N 


9.067 


31.1 


N 


488 


3.338 ± 1.850 


0.000 ± 0.000 


3.338 





8.614 


2.8 


N 


230 


0.000 ± 0.000 


0.263 ±0.169 


0.263 


N 


7.620 


27.0 


N 


409 


0.002 ± 0.001 


0.002 ± 0.001 


0.000 





7.574 


2.2 


N 


94 


1.572 ± 1.097 


7.631 ± 2.488 


6.058 


* 


7.466 


28.5 


N 


344 


2.701 ± 0.497 


2.088 ± 0.515 


0.613 




7.465 


15.3 


N 


130 


0.099 ± 0.047 


0.221 ± 0.069 


0.122 


* 


7.054 


31.7 


N 


111 


1.002 ± 0.323 


0.000 ± 0.000 


1.002 





6.985 


11.4 


N 


109 


0.495 ± 0.322 


1.318 ±0.852 


0.823 




6.865 


18.7 


N 


42 


1.144 ±0.362 


0.083 ± 0.389 


1.061 


* 


6.829 


0.6 


N 


63 


0.000 ± 0.000 


0.424 ±0.143 


0.424 


N 


6.451 


17.4 


N 


12 


2.297 ±0.319 


1.505 ±0.331 


0.792 


* 


6.159 


21.7 


N 


469 


0.088 ± 0.073 


0.000 ± 0.000 


0.088 





6.107 


18.1 


N 


144 


0.176 ± 0.297 


0.751 ±0.361 


0.575 




6.087 


22.8 


N 



Table 7: Same as Table 5 



IAU designation 


INPOPlOc 


INPOPlOd 


Diff 


S 


Impact 


L 




number 


10 12 xMq 


10 12 xMq 


10 12 xM G 




m 


% 




356 


4.173 ± 0.868 


4.173 ±0.902 


0.000 




5.759 


2.2 


N 


712 


0.000 ± 0.000 


1.228 ±0.267 


1.228 


N 


5.745 


2.2 


N 


88 


1.340 ±0.866 


0.000 ± 0.000 


1.340 





5.742 


1.4 


N 


60 


0.402 ± 0.221 


0.282 ± 0.268 


0.120 




5.733 


3.8 


N 


50 


0.686 ±0.187 


1.031 ±0.566 


0.345 




5.702 


2.0 


N 


128 


4.699 ± 1.522 


0.000 ± 0.000 


4.677 





5.624 


3.0 


N 


5 


0.448 ±0.165 


0.913 ± 0.220 


0.466 


* 


5.533 


15.8 


N 


59 


4.332 ± 0.607 


1.364 ± 1.097 


2.968 


* 


5.325 


12.1 


N 


98 


1.414 ± 0.603 


2.100 ±0.705 


0.686 




5.195 


15.1 


N 


194 


6.387 ± 0.701 


4.380 ±0.819 


2.007 


* 


5.145 


2.9 


N 


51 


3.546 ± 0.748 


3.639 ± 0.937 


0.093 




5.109 


15.6 


N 


156 


3.263 ± 0.438 


3.089 ± 0.576 


0.174 




5.103 


19.3 


N 
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Table 8: Asteroid masses found in the recent literature compared with the values estimated in INPOPlOc and INPOPlOd. The 
uncertainties are given at 1 published sigma. 



IAU designation 


INPOPlOc 


Close-encounters 


Refs 


INPOPlOd 




number 


10 12 


xM Q 


10 12 x M Q 




10 12 xM Q 


% 


5 


0.448 


±0.165 


1.705 ±0.348 


Zielenbach (2011) 


0.913 ± 0.220 


15.8 


12 


2.297 


±0.319 


2.256 ± 1.910 


Zielenbach (2011) 


1.505 ±0.331 


21.7 


20 


0.000 


± 0.000 


1.680 ±0.350 


Baeretal. (2011) 


1.921 ±0.661 


39.9 


24 


7.6 


± 1.6 


2.639 ± 1.117 


Zielenbach (2011) 


2.2 ± 1.7 


36.0 


27 


0.000 


± 0.000 


1.104 ±0.732 


Zielenbach (2011) 


1.511 ±0.982 


29.5 


45 


3.523 


±0.819 


2.860 ± 0.060 


Marchis et al. (2008) 


1.518 ±0.962 


21.0 


59 


4.332 


± 0.607 


1.448 ±0.0187 


Zielenbach (2011) 


1.364 ± 1.097 


12.1 


94 


1.572 


± 1.097 


7.878 ±4.016 


Zielenbach (2011) 


7.631 ±2.488 


28.5 


130 


0.099 


± 0.047 


3.320 ± 0.200 


Marchis et al. (2008) 


0.221 ± 0.069 


31.7 


139 


0.000 


± 0.000 


3.953 ± 2.429 


Zielenbach (2011) 


3.579 ± 0.595 


32.0 



20 
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Appendix A: Analytical solution 

This section presents the analytical solutions of Equation 1. 
Let, I\ and h be the integral solutions of Equation 1 using 
Equation 3 and 4 i.e. , 



From the geometry of Figure 3, the lower limit of Equation 
A. 3 can be written as 



arctan 



-R S /E cos a 



arctan 



-Rs/e cos a 



R 



s/e sin a 



h 



- f * af-V 



and 



h 



r Kf 



dL 



+ Bl±- 



-2 



dL . 



From the geometry (Figure 3) we define 
P — Rs/e sin a , 
Ldc = Re/sc ~ Rs/e cos a , 



^DE - 1 r 



I 2 = L 2 +R 2 S/E - ILRs/E cos a 



(A.l) 

with 

(A. 2) cot a = tan 
Hence, 
arctan 



arctanl - cot a 



Rs/e cos a\ n 
— a 



Similarly, the upper limit of Equation A. 3 can be written as 



arctan 



with 



Reisc - Rs/e cos a 



— arctan 



Re/sc - Rs/e cos a 
Rs /e sin a 



where a and ft are the angle between the Sun-Earth-Probe 
(SEP) and the Earth-Sun-Probe (ESP). P is the MDLOS from ^ _ Rs/£ CQS Q = Rs/£ ^ ff 
the Sun. With these expressions, I\ can be written as 



tan f} 



~-BR% f 

JL, 



BR" 



L, lc 



dL 



(L 2 + R 2 S/E -2L R s/E cos ayl 2 



Hence, 



arctan 



dL 



([L - R S /e cos a] 2 + R 2 SIE sin 2 ayl 2 



Re/sc - Rs/e cos a 



Now Equation A. 3 is given by 



Assuming, 
x — L — Rs/e cos a , 
a = Rs/e sin a , 
dx = dL , 

with L — at the Earth station (L£ art /, l/n ) and L = Re/sc at 
the spacecraft (L s / C ). Then the integral I\ can be written as 



BRl 

a e 

with 



•/3+a-$ 



(tan 2 6+iy/ 2 



dO 



sec 6 = tan 6+1 
and 

cos 6 ■■ 



1 



h =BR 



Rs/e cos a 

L 



Rsie cos a 



dx 



sec 6 

Therefore, the integral I\ can be written as 



B Rq r s £/sc-Ss/£ cos a 



@ J -Rs/e cos a 

Now let 

x 

— = tan 6 , 
a 

and 

dx — a sec 2 6 d6 . 
Therefore, 
BR 



(x 2 + a 2 y' 2 

dx 

(i + 4) e/2 ' 



br% r p+ - 

Ja-f 



(cos ey- 2 de 



(A.4) 



The maximum contribution of the integral occurs at = 0. 
To solve this integral, Taylor series expansion was used and for 
9 near zero, it can be given as 



f(0)=f(O) + 



df 
dO 



2 (d 2 f\ 3 (d 3 f 



2\\d0 2 3! \ dO 3 



• + 9( „ 



4! [de* 



(A.5) 



' i 



5^1 r 

® J an 



(R E/SC- R S/E cos a) 



a sec 2 
(tan 2 0+iy' 2 



with 



dO 



(A3) 



f(0) = (cos 0) 



-2 



Then 
/(O) = 1 

^ = -(e-2)(sin0)(cosO) { ^ ) 

m -o 
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3e 2 - 14e+ 16 



120 



08 + or - n/2Y -(a- tt/2)- 



21 
(A.6) 



Now, Equation A. 2 can be written as h- ha + hb where, 



ha = f * a(-L) dL , 



<d 2 f 



dff- 



5- = (e - 2) (e - 3) (sin 1 0) (cos 9) 



(er-4) 



and 



hb 



f* 44-) 



dL 



(A.7) 



(A.8) 



-(e - 2) (cos 0) (e - 2) 
= -(e-2) , 



1=0 



iff 

[de 2 



Using a similar approach to the previous integral, one can 
write 



ha 



ARt rP +a 



AR 4 & 



2 9d0 



+(e - 2) (3e - 8) (sin 0) (cos 0) 



(e-3) 



r - 2 

cos^ 

J a- % 



-(2 0+ sin 29) 
4 





(d 4 f 



ARi 



0=0 



4 a 3 



(2/3- sin 2(a + B) + sin 2a) 



,d9* 



(e - 2) (e - 3) (e - 4) (e - 5) (sin 4 0) (cos 0) 



(e-6) 



ARi 



4 Rl/ E sin?a 



(2/3- sin 2(a + /3) + sin 2a) 



-(e - 2) (e - 3) (6e - 20) (sin 2 0) (cos 0) 



(E-4) 



Similarly, Equation A.8 can be written as 



<d*f 

yd& 



+(e-2)(3e-8)(cos6) 



(e-2) (3e-8) 



(e-2) 



hb 



Now, Equation A. 5 can be written as 



(cos 9) {e - 2) = 1 - (e - 2) — + (e - 2) (3e - 8) — + .... 9 (0") 



BRl 



rP+a 



dO 



2! 



1 



e-2 



3e 2 - 14e+ 16 



4! 



e 4 + ....9(9") 



BRl 



Rs ie sin a 



2 24 
By neglecting the higher order terms, the integral (Equation 



Now by using ha and hb, expression h can be written as 



A.4) can be written as 



/1 



h = 



AR« 



BR% j 

a- 1 J 


-i (' 


BR% 


e-2 


a- 1 


r 6 


BR% 


e-2 

r 6 


36 2 - 


14e+ 16 / 



9 l + 



24 



4 R\i E sin 3 a 



(2/3 - sin 2(a + B) + sin 2a) 



B Ri 
+ —B . 

Rs/e sin a 



dO 



(A.9) 



9' + 



3e 2 - 14e+ 16 



120 



5 



120 



(J3 + a - tt/2) 3 - (a - n/2)- 



(B + a- tt/2) 5 -(a- tt/2) 5 



By substituting the value of a, we can write the integral I\ 



as 

h = 



BR" 



(R S /e sina) 



E-l 



B- e —^((J3 + a- tt/2) 3 -(a- ;r/2) 3 ) 



